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Abstract 

To investigate the cosmic no hair conjecture, we analyze numerically 1-dimensional 
plane symmetrical inhomogeneities due to gravitational waves in vacuum spacetimes with 
a positive cosmological constant. Assuming periodic gravitational pulse waves initially, 
we study the time evolution of those waves and the nature of their collisions. As mea- 
sures of inhomogeneity on each hypersurface, we use the 3-dimensional Riemann invariant 
X = ( 3 'Rijki ^K l i kl and the electric and magnetic parts of the Weyl tensor. We find a tem- 
poral growth of the curvature in the waves' collision region, but the overall expansion of the 
universe later overcomes this effect. No singularity appears and the result is a "no hair" de 
Sitter spacetime. The waves we study have amplitudes between 0.020A < X 1 / 2 < 125. OA 
and widths between 0.080Z# < I < 2.5Z#, where lu = (A/3)" 1 / 2 , the horizon scale of de 
Sitter spacetime. This supports the cosmic no hair conjecture. 
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The widely accepted inflationary universe model is an attempt to solve many long- 
standing cosmological difficulties — such as the flatness, horizon, and monopole problems — in 
the standard big-bang model. The explanation of the present isotropy and homogeneity of 
the universe are one of those problems. The inflationary scenario may resolve this problem by 
introducing a rapid cosmological expansion. But the scenario is usually studied in an isotropic 
and homogeneous spacetime, so that to determine whether inflation really isotropizes and ho- 
mogenizes the universe and to see whether inflation is generic, we have to study an inflationary 
model in anisotropic and inhomogeneous spacetimes as well. 

To explain the generality of inflation, the so-called cosmic no hair conjecture has proposed |2|]. 
The conjecture is: "All initially expanding universes with positive cosmological constant A ap- 
proach the de Sitter spacetime asymptotically". There exist, however, simple counter examples. 
For example, a closed Friedmann-Robertson- Walker spacetime with a positive A recollapses if 
the initial energy density is large enough. On the other hand, the de Sitter spacetime is stable 
against perturbations and many other models also support this conjecture, so that we expect 
we can prove a suitably refined version of this conjecture. 

For homogeneous but anisotropic spacetimes, Wald|| showed that, except in the Bianchi 
IX case, initial expansion and positive A force an approach to the de Sitter spacetime within 
one Hubble expansion time th = ^/A/3. And for Bianchi IX, the additional initial condition 
that A be larger than one third of the maximum value of 3-dimensional Ricci scalar leads to 
the same conclusion ||, ||. For inhomogeneous spacetimes, on the other hand, the only practical 
method we have at present is to solve the Einstein equations numerically. Assuming a spher- 
ically symmetric spacetime, Goldwirth and Piran|| studied the behavior of inhomogeneous 
distributions of a scalar field. They found a sufficient criterion for inflation that the physical 
scale of inhomogeneity of the scalar field be larger than several times the horizon scale. 

We may, however, wonder whether this result is valid for cases with other symmetries or 
for more generic spacetimes. Assuming a plane symmetric spacetime, Kurki-Suonio et. a/.|| 
studied new inflationary model and found that inflation occurs only when the potential is flat 
enough. The spacetime they studied is, however, not sufficiently inhomogeneous in order to 
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investigate the "no hair conjecture" . This is the first reason why we study a plane symmetric 
spacetime in this paper. 

Besides the inflaton field, there exist other sources of inhomogeneity, one of which is gravita- 
tional waves. Although any linear gravitational waves in de Sitter space will decay exponentially 
and desappear soonQ, we do not know what will happen if strong gravitational waves exist. If 
the strong gravitational waves are localized, those may collapse into black holes. Even if those 
are not localized as a plane symmetric case, such a strong gravitational field may evolve into a 
naked singularity. In fact, Szekeres and Khan- Penrose |J showed that a collision of two plane 
symmetric gravitational waves in Minkowski background spacetime will form a singularity. We 
may wonder what will happen with strong gravitational waves in de Sitter background. A sin- 
gularity may be formed and prevent inflation. In this report, we study such an inhomogeneity 
due to gravitational waves. Because they cannot exist in a spherically symmetric spacetime, 
gravitational wave inhomogeneities may play a different role in the homogenization process 
from that of the inflaton field analyzed previously. This is the second reason why we study a 
plane symmetric spacetime here. 

Recently, Nakao et al. examined time symmetrical initial data for Brill waves in an axial 
symmetrical spacetime with cosmological constant and found that waves with large gravita- 
tional mass does not always provide trapped surfaces citenakaogw. They also found that a dust 
sphere with large gravitational mass in a background de Sitter spacetime does not collapse 



to form a black hole spacetime [10], and that there exists an upper bound on the area of the 



apparent horizon of a black hole in an asymptotically de Sitter spacetime [H], so that we may 
conjecture that a large inhomogeneity does not necessarily prevent cosmological constant-driven 
inflation. From these results, it seems plausible that in an inflationary era, inhomogeneities 
will simply evolve into many small black holes in a background de Sitter universe; the only 
worrisome possibility is that a naked singularity might form. Disposing of this worrisome pos- 
sibility may thus become the main problem in the study of the cosmic no hair conjecture. Our 
1-dimensional case seems well suited to address this problem. This is the third reason since a 
singularity formed in a plane symmetric spacetime is always naked. 
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For the above three reasons, we have carried out a numerical study of inhomogeneities due 
to plane-symmetric gravitational waves with a cosmological constant. Since we treat a plane 
symmetric spacetime and gravitational waves, our simulation may be thought of a complemen- 
tary work to the spherically symmetric one discussed by Goldwirth and Piran. Since we have 
not studied an inflaton field in this paper, however, our result may not be applied to a realistic 
universe as it is. (See also our comment in the references). 

To analyze the behavior of 1-dimensional inhomogeneities due to gravitational waves, we 
assume spacetime is plane symmetric and vacuum. We also assume that a positive cosmological 
constant A exists. Under these assumptions, we examine whether such a spacetime leads to an 
inflationary era, and whether such initial inhomogeneities and anisotropies smooth out during 
inflation periods. 

We use the ADM formalism to solve the Einstein equation 

Rfiu — -^9tivR + A# Mi/ = 0, labeleinstein (1) 

with the metric 

ds 2 = -(a 2 - (3 2 h xx )dt 2 + 2f3dtdx + ^ xx dx 2 + ^ yy dy 2 + ^ yz dydz + ^ zz dz 2 : , (2) 

where the lapse function a, the shift vector $ = (/3, 0, 0), and the 3-metric 7^ depend only on 
the time t and propagation direction coordinate x. 

We use the Hubble expansion time th = ^A/3 as our time unit, which is a characteristic 
expansion time of the expanding Universe. And our unit of length is also normalized to the 
horizon length of the de Sitter universe Ih = 1/ ^A/3. 

Our simulation procedure follows Nakamura et al. ||12|| : (i) Determine initial values by solving 
two constraint equations, (ii) Evolve time slices by using the dynamical equations, (iii) Check 
the results of (i) and (ii) using the two constraint equations on every time slice. 

To determine the initial data, we use the 5th order Runge-Kutta method (Fehlberg method), 
and to integrate dynamical equations, we use a finite difference method with 400 meshes. In 
our calculation, the maximum errors in the Hamiltonian and momentum constraint equations 
are O(10~ 4 ) and O(10~ 6 ) on the initial hypersurface, respectively, and these accuracies are 
maintained even after evolution. 
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We use the York-O'Murchadha's conformal approach [[13] to get initial values. Setting 7^ = 
4> 4/ yij, the Hamiltonian constraint equation becomes 

8 (3) A0 = {3) R<f) - AjA^^- 7 + 2{-K 2 - A)<f) 5 (3) 

3 

where is the 3-dimensional Laplacian of a conformal metric 7^, is the 3-dimensional 
scalar curvature, hatA^ is the trace-free part of the extrinsic curvature and hatK(= K) is 
a trace part of K^. The quantities with a hat denote physical variables in the conformal frame. 
In this approach, jij,K and the transverse- traceless(TT) part of the extrinsic curvature Kf 3 T 
are left to our choice. 

We assume constant mean curvature on the initial hypersurface, 

K = -v / 3A(l + 5 K ) = const. on Sigma(t = 0), (4) 

where 5 k is a constant. Here, 5k is introduced to find the periodic gravitational waves. 5k = 
is inconsistent with the presence of gravitational waves in the initial data. The reason is, when 
we have pulse-like inhomogeneities under periodic boundary condition, there are two energy 
sources; one is a positive cosmological constant and the other is the energy of pulse waves. 
These energy give contributions to the expansion of the universe, i.e. the first term in the r.h.s. 
of (|4]) is due to the cosmological constant, while the second one (8k) is due to the pulse waves. 

With this slicing, the momentum constraint equation becomes trivial and A^ = K^ T . If 
we also impose — > 1 at the numerical boundary edges, the value of 5k is fixed; it is about 
O(10~ 2 ) in the physical frame. 

For coordinate conditions, we impose a geodesic slicing condition such as a = 1 and (3 = 0. 
This slicing condition entails the risk that our numerical hypersurface may hit a singularity and 
stop there. If no singularity appears, however, then it may be the best coordinate condition for 
revealing whether de Sitter space emerges as a result of time evolution. 
As for the initial data, we consider two complementary cases: 
[case 1] The inhomogeneities reside in the 3-metric garhma i j and Kj^ = 0, 
[case 2] The 3-metric is conformally flat and the inhomogeneities reside in Kf? ' . 
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In this Brief Reports, we present mainly the results of simulations for the [case 1], and we 
examine a pulse-like wave propagating in the x-direction, and expressed by the metric in York's 
conformal frame as 

diag(7, J ) = (l,l + ae-^ 2 ,l), (5) 

where a and x are free parameters. Varying these parameters, which we interpret as the 
wave's amplitude and width, we study the effects of a family of gravitational pulse waves on 
the spacetime structure. 

As a measure of inhomogeneity, we first use the Riemann invariant scalar X = ^Rijki ^R^ kl , 
where ^Rijki is the Riemann tensor of the 3-metric on E, and we introduce a dimensionless 
variable normalized by the cosmological constant, 

jl/2 

C{t,x) = — onS(t), (6) 

and call it the "curvature" hereafter. We estimate the magnitude of the inhomogeneities 
in the 3-space E by the maximum value of this "curvature" on each slice, i.e. C max (t) = 
max{C(t,x) | zgs}. 

As we mentioned a moment ago, a pulse-like wave has two characteristic physical dimensions, 
a width and an amplitude. We use C max (t) as an amplitude measure, and we define the width 
l(t) by the proper distance between two points where the trace 7 of 3-metric 7^ (the square 
of the 3- volume) decreases by half from its maximum value ^max- When a and xq are given, 
the initial form of a gravitational wave is determined, and the width 1(0) and C max (0) may be 
calculated. 

In Figure 1, we show a typical time evolution of this model for the case a = — 0.10, xo = 
0.20; I = 0.33, C max (0) = 1.25. We see that the wave given on the initial surface propagates both 
in the +x and — x directions, and the "curvature" seems to be superposed in the interacting 
region of the waves. The spacetime, however, finally succumbs to the overall expansion driven 
by the cosmological constant, and becomes virtually indistinguishable from de Sitter spacetime 
within one Hubble expansion time r#. 

We also use 4-dimensional variables as a measure to see the homogenization process. In 
the vacuum spacetime, if the Weyl tensor r h oa vanish, the spacetime is homogenized and 
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isotropized. We use the decomposition of the Weyl tensor: E pa = C ppau n ,1 n L ' ', B pa = *C ppav n^n v ', 
where *C pupa = \C a ^ p sigma £ a /3 P u and is a timelike vector orthogonal to the hypersurface E. 
In analogy to electromagnetism, the 3-dimensional variables E pa and B pcr are called an electric 
and a magnetic parts of the Weyl tensor, and we can reconstruct the Weyl tensor completely 
from this pair of tensors. We introduce 

H(t,x) = E pu E^ + B pu B^, (7) 

as a sort of gravitational "super-energy" (in fact the purely timelike component of the Bel- 
Robinson tensor PHI). If TC(t,x) — > 0, then C^ upa — > unless the hypersurfaces become null. 
We find that H(t, x) almost vanishes within one Hubble expansion time, and that this indicates 
that our spacetime does indeed approach a homogeneous one (i.e. de Sitter spacetime). 

In our simulation, I and C max (0) range between 0.080£# < I < 2.5Z# and 0.020 < C max (0) Zeg80.0. 
From the present results, we conclude that for any large Riemann invariant and/or small width 
inhomogeneity on the initial hypersurface, the nonlinearity of the gravity has little effect and 
the spacetime always evolves into a de Sitter spacetime. 

Jensen and Stein-Schabes[|15[ proved an extension of Wald's theorem, showing that inflation 
invariably eliminates all inhomogeneities consistent with a non-positive scalar 3-curvature. Such 
inhomogeneities are hardly generic, however. In our models, we find that is initially 

negative, but it eventually becomes positive along the propagations of the distortions and 
finally ( 3 >R — > when the spacetime is homogenized (see Figure 2). This convinces us that 
the behavior we have found is not the case envisioned by Jensen and Stein-Schabes. 

Next, we simulated collisions of two gravitational pulse waves. The motivation for such a 
simulation is the existence of exact solutions for colliding plane waves in which scalar curvature 
singularities occur. Szekeres and Khan- Penrose || were the first to find such solutions. Szekeres 
treats the collision of two thick or sandwich plane waves, while Khan-Penrose treat two 5 
function-like shock waves. These solutions show the nonlinearity of the gravity. Thus we 
interest in such a phenomena in more general and realistic situations. 

It may be interesting to see whether singularities also occur when gravitational waves collide 
in the presence of a cosmological constant. To investigate this possibility, we construct initial 
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data representing two nearby pulse waves. Since we impose periodic boundary conditions, there 
is no difference between the present case and the case with a single pulse waves in one half of 
our numerical region. But in the present setting the collision does not occur on the boundary, 
so it may be better to see the interaction of two waves. In Figure 3, a typical example of 
the time evolution of the "curvature" is shown. We see that the waves start to propagate in 
opposite directions (+x and —x) and then collide with each other. When the waves collide, 
the "curvature" is just superposed as in Fig.l, but more clearly. After the wave-crossing, each 
wave propagates away. (This behavior is reminiscent of solitonic waves.) Consequently, the 
spacetime is again homogenized by the expansion of the universe within one Hubble expansion 
time tauH- For a wide range of initial widths and "curvatures" (0.080Z# < I < O.IOIh, 
40.0 < C max (0) < 125.0 and the periodic distance d is 0.20l H < d < 0.501 h), all inhomogeneities 
decay below 1 % of their initial "curvatures" within one Hubble expansion time. 

That these spacetime become homogeneous even if we include colliding waves is consistent 
with the calculations by Centrella and Matzner[[ITJ. They examined the collision of gravita- 
tional shock waves in an expanding Kasner background both analytically and numerically, and 
concluded that such a collision leads to no singularity. The expansion due to the cosmological 
constant in our simulation is exponential — much faster than Kasner (power-law) expansion. 
Thus, our result would be expected from their result, although here we investigated pulse 
wave collisions systematically for a wide range of initial parameters and extracted the generic 
behavior of gravitational waves in expanding universe. 

We also computed the evolution of [case 2] initial data, 7^ = 5^ and diag ( hatKj^) = 
(0, a q~( x / x o) 2 ^ — a e -( x / x o) 2 ^ anc [ found that all inhomogeneities present in them also decay 
within one Hubble expansion time. We also tried similar calculations with a different shape 
for the initial pulse waves, i.e., assuming a [1 + acos 2 (x/xo)]-form for the ^ yy components, and 
found the same results for time evolution. We conclude from this that our result is generic: all 
initial inhomogeneities decay within one Hubble expansion time and disappear. 

The present result suggests no additional condition of the kind needed to make cosmic no 
hair conjecture a theorem, but it does give a wide class of examples in which gravitational 
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inhomogeneities succumb to cosmological constant- driven expansion. We could not find any 
characteristic scale length. All inhomogeneous spacetimes expand to de Sitter universe. In this 
sence, our results are different from those of Goldwirth and Piran||, who simulated effects of 
spherically symmetric inhomogeneities due to a realistic inflaton field. We do not yet know 
whether our results reflect the generic behavior of inhomogeneities driven by gravitational 
waves, or whether they are largely due to our ansatz of plane symmetry. In order to answer 
these questions, further simulations including a inflaton field, or more gravitational degrees of 
freedom — for example, those present in cylindrical or axially symmetrical, or even more general 
spacetimes — -are required. Such simulations are underway[T7| . 
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Figure Captions 

Figure 1: A typical example of the time evolutions of propagating plane waves. The "cur- 
vature" C(t,x) (see eq. (|)) is shown for the initial pulse shape of I = 0.33/# 
and C max (t = 0) = 1.25 (a = -0.10, x = 0.20 in eq. (§)). We see the initial 
pulse propagates both in the +x and —x directions, but is smoothed out by the 
expansion of the universe within one Hubble expansion time. 

Figure 2: The evolution of ( 3 'R at x = (center) for the same initial data as in Fig.l. This 
shows also the maximum deviations from zero at each time. We see that ®R 
does not remain negative for all time, so that our simulations are not governed 
by the case proven by Jensen and Stein-Schabes. 

Figure 3: The time evolution of the "curvature" C(t, x) resulting from waves that are located 
closely. Two waves are the same form, I = 0.10Z# and C max (t = 0) = 51.0, and the 
periodic distance is 0.30Z#. We see them collide and in the collision region, the 
"curvature" seems to be superposed, but finally the spacetime is homogenized by 
the expansion of the universe. 
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